Search for: All records

The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (FIR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine-structure lines and emission features from dust and ionised and neutral atomic and warm molecular gas in the rest-frame mid-infrared (MIR) and FIR provide unmatched diagnostic opportunities to determine the properties of gas and dust, measure gas-phase metallicities, and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low- and high-resolution FIR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a 1.8 m diameter, cryogenically cooled FIR observatory currently at the conception stage, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZand the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to 10¹¹ L_⊙beyond cosmic noon and at least up toz = 4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbons (PAHs), and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow us to trace both the star formation and black hole accretion rates (SFRs and BHARs), the gas-phase metallicities, and the mass-outflow rates of cold gas in hundreds to thousands of individual galaxies toz = 2. 

We present a study of new 7.7–11.3 μm data obtained with theJames WebbSpace Telescope Mid-InfraRed Instrument in the starburst galaxy M 82. In particular, we focus on the dependency of the integrated CO(1–0) line intensity on the MIRI-F770W and MIRI-F1130W filter intensities to investigate the correlation between H₂content and the 7.7 and 11.3 μm features from polycyclic aromatic hydrocarbons (PAH) in M 82’s outflows. To perform our analysis, we identify CO clouds using the archival¹²CO(J = 1 − 0) NOEMA moment 0 map within 2 kpc from the center of M 82, with sizes ranging between ∼21 and 270 pc; then, we compute the CO-to-PAH relations for the 306 validated CO clouds. On average, the power-law slopes for the two relations in M 82 are lower than what is seen in local main-sequence spirals. In addition, there is a moderate correlation betweenI_{CO(1 − 0)} − I_7.7 μm/I_11.3 μmfor some of the CO cloud groups analyzed in this work. Our results suggest that the extreme conditions in M 82 translate into CO not tracing the full budget of molecular gas in smaller clouds, perhaps as a consequence of photoionization and/or emission suppression of CO molecules due to hard radiation fields from the central starburst. 

We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain-size distributions, and ionization states with a new model for dust evolution in galaxy simulations. We apply these models to three idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and a starburst disk. Our main results are as follows. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to shatter grains efficiently, driving up the fraction of ultrasmall grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain-size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to investigate fully the complex interplay of processes that drive PAH band luminosities in galaxies. 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

Abstract Accurately determining gas-phase metal abundances within galaxies is critical as metals strongly affect the physics of the interstellar medium. To date, the vast majority of widely used gas-phase abundance indicators rely on emission from bright optical lines, whose emissivities are highly sensitive to the electron temperature. Alternatively, direct-abundance methods exist that measure the temperature of the emitting gas directly, though these methods usually require challenging observations of highly excited auroral lines. Low-lying far-infrared (FIR) fine structure lines are largely insensitive to electron temperature and thus provide an attractive alternative to optically derived abundances. Here, we introduce the far-infrared abundance (FIRA) project, which employs these FIR transitions, together with both radio free–free emission and hydrogen recombination lines, to derive direct, absolute gas-phase oxygen abundances. Our first target is M101, a nearby spiral galaxy with a relatively steep abundance gradient. Our results are consistent with the O ++ electron temperatures and absolute oxygen abundances derived using optical direct-abundance methods by the CHemical Abundance Of Spirals (CHAOS) program, with a small difference (∼1.5 σ ) in the radial abundance gradients derived by the FIR/free–free-normalized versus CHAOS/direct-abundance techniques. This initial result demonstrates the validity of the FIRA methodology—with the promise of determining absolute metal abundances within dusty star-forming galaxies, both locally and at high redshift. 

Abstract We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral field spectroscopy of the nearby luminous infrared galaxy NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST/NIRSpec to study the 3.3μm neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ∼200 pc scales. A clear change in the average grain properties between the star-forming ring and the central AGN is found. Regions in the vicinity of the AGN, with [Neiii]/[Neii] > 0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios, indicating that smaller grains are preferentially removed by photodestruction in the vicinity of the AGN. PAH emission at the nucleus is weak and shows a low 11.3/3.3 PAH ratio. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3μm–to–total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on subkiloparsec scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3μm PAH-derived star formation rate (SFR) in the ring is 27% higher than that inferred from the [Neii] and [Neiii] emission lines, the integrated SFR derived from the 3.3μm feature would be underestimated by a factor of 2 due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high redshift. 

Abstract On 11 September 2021, two small thunderstorms developed over the Telescope Array Surface Detector (TASD) that produced an unprecedented number of six downward terrestrial gamma ray flashes (TGFs) within one‐hour timeframe. The TGFs occurred during the initial stage of negative cloud‐to‐ground flashes whose return strokes had increasingly large peak currents up to 223 kA, 147 GeV energy deposit in up to 25 1.2 km‐spaced surface detectors, and intermittent bursts of gamma‐rays with total durations up to 717 s. The analyses are based on observations recorded by the TASD network, complemented by data from a 3D lightning mapping array, broadband VHF interferometer, fast electric field change sensor, high‐speed video camera, and the National Lightning Detection Network. The TGFs of the final two flashes had gamma fluences of and 8, logarithmically bridging the gap between previous TASD and satellite‐based detections. The observations further emphasize the similarity between upward and downward TGF varieties, suggesting a common mechanism for their production. 

Abstract Optical emissions associated with Terrestrial Gamma ray Flashes (TGFs) have recently become important subjects in space‐based and ground‐based observations as they can help us understand how TGFs are produced during thunderstorms. In this paper, we present the first time‐resolved leader spectra of the optical component associated with a downward TGF. The TGF was observed by the Telescope Array Surface Detector (TASD) simultaneously with other lightning detectors, including a Lightning Mapping Array (LMA), an INTerFerometer (INTF), a Fast Antenna (FA), and a spectroscopic system. The spectroscopic system recorded leader spectra at 29,900 frames per second (33.44 s time resolution), covering a spectral range from 400 to 900 nm, with 2.1 nm per pixel. The recordings of the leader spectra began 11.7 ms before the kA return stroke and at a height of 2.37 km above the ground. These spectra reveal that optical emissions of singly ionized nitrogen and oxygen occur between 167 s before and 267 s after the TGF detection, while optical emissions of neutrals (H I, 656 nm; N I, 744 nm, and O I, 777 nm) occur right at the moment of the detection. The time‐dependent spectra reveal differences in the optical emissions of lightning leaders with and without downward TGFs.